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PREFACE 

The  work  described  in  chls  publication  was  performed  by  the 
Mathematics  Clinic,  Claremont  McKenna  College,  under  contract  to  the 
Jet  Propulsion  Laboratory,  an  operating  division  of  the  California 
Institute  of  Technology.  This  activity  Is  sponsored  by  the  Jet 
Propulsion  Laboratory  under  contract  NAS7-918,  RE182,  A187  with  the 
National  Aeronautics  and  Space  Administration,  for  the  United  States 
Army  Intelligence  Center  and  School. 
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I  INTRODUCTION  AND  OVERVIEW 


This  is  the  final  report  of  the  1935-36  Claremont  Graduate  School  and 
Claremont  McKenna  College  Mathematics  Clinic  concerning  a  problem  proposed 
by  the  Intelligence  Analysis  Group  at  the  Jet  Propulsion  Laboratory.  The 
problem  is  to  determine  confidence  regions  for  the  location  of  objects  or 
emitters  on  which  bearings  are  taken  from  two  or  more  sensors  whose 
positions  are  known.  -- 

N 

This  problem  for  the  two-sensor  case  was  studied  by  a  different  clinic 
during  the  1  934  -85  academic  year.^This  clinic  studied  the  classical 
approach  to  determining  such  regions  which  involved  the  assumption  that  an 
error  of  observation  of  the  line  of  bearing  displaces  the  line  parallel  to 
itself.  By  dropping  this  assumption,  they  found  that  the  classical 
probability  regions,  which  are  ellipses,  conta  ined.  on!  y  50*  to  SO*  of  the 
probability  claimed.  One  goal  of  this  year’s  clinic  was  to  determine  the 
reason  for  this  very  large  discrepancy.  As  it  turned  out,  the  reason  was 
neither  the  dropping  of  the  parallel  displacement  assumption,  nor 
computational  errors  in  the  clinic's  work,  but  two  equations  in  the  original 
publication  describing  the  classical  approach  were  printed  incorrectly. 
These  errors  are  probably  transcription  errors,  but  no  corrections  or 
references  to  them  were  found  in  a  search  of  the  relevant  literature.  ^ 

In  Section  II,  the  two-sensor  problem  is  described.  The  parallel 
displacement  assumption  is  dropped.  Based  on  tne  assumption  that  errors  of 
the  reported  angles  are  normally  distributed,  tne  probability  density 
function  of  the  points  of  intersection  of  the  computed  lines  of  bearing  is 
given. 

In  Section  III  and  the  Appendix,  a  re-der  i  vat  ion  of  the  equations  of 
the  classical  model  is  carried  out  and  the  correct  equations  for  the 
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probability  regions  (called  Stansfield  ellipses)  are  determined. 

In  Section  IV,  it  is  shown  that  the  Stansfield  ellipses  as  corrected  by 
the  clinic,  actually  contain  the  claimed  probabilities  with  only  very  minor 
differences.  However,  the  densities  above  and  below  the  location  of  the 
emitter  are  different  due  to  the  dropping  of  the  parallel  displacement 
assumption.  The  contour  maps  shown  in  figure  5  show  this  clearly. 

In  Sections  V  and  VI,  two  new  probability  regions  are  developed,  the 
"Clinic  region"  and  the  "Image  region".  These  are  attempts  at  determining 
probability  regions  which  contain  the  same  probability  as  a  Stansfield 
ellipse,  but  whose  areas  are  smaller  than  the  area  of  the  ellipse. 

In  Section  VII,  the  Clinic,  Image  and  Stansfield  regions  are  analyzed. 
It  is  shown  that  the  Image  region  is  not  as  "efficient"  (larger  area  for 
given  probability)  as  the  Clinic  and  Stansfield  regions.  The  Clinic  region, 
on  the  other  hand,  is  smaller  than  the  Stansfield  ellipse  for  a  given  level 
of  probability,  but  the  difference  is  so  small  that  it  is  hardly  worth  the 
effort  to  determine  this  much  more  complicated  region.  However,  the 
geometric  center  of  the  Clinic  region  moves  from  below  to  above  the  true 
location  of  the  emitter  as  the  desired  probability  of  the  region  is 
increased  (see  figure  13),  wnile  the  geometric  center  of  the  Stansfield 
ellipse  remains  fixed.  This  result  may  be  of  practical  value. 

In  Section  VIII,  the  clinic  approaches  the  three-sensor  problem  in  a 
manner  similar  to  the  two-sensor  analysis.  The  purpose  is  to  find  the  true 
density  function  when  the  parallel  displacement  assumption  is  dropped. 
Several  approaches  are  taken  to  obtain  a  discrete  probability  density 
function  to  approximate  the  true  density  function,  which  is  continuous.  The 
different  approaches  produce  very  different  results.  More  .._rk  should  be 
done  in  this  area. 
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Section  IX  is  very  different  from  the  other  sections.  Here,  a  model 


for  a  single,  active  sensor  problem  is  developed.  An  active  sensor  measures 


both  direction  and  distance,  taking  readings  from  two  different  locations. 


Under  the  assumption  that  the  systematic  errors  (direction  and  distance)  are 


the  only  significant  errors,  a  scheme  to  determine  the  true  location  of  the 


emitter  is  presented. 


and  Og*.  He  assume  the  two  sensors  ai 
that  tne  pair  of  readings  has  the  biv; 
probability  density  function  for  (A,U 


Ho, 3)  = 
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exp 
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where  the  values  of  the  random  variab 


<*0  and  B0  are  the  true  angles  of  bear 
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0  *  tan ( 3 ) 


x  = 


tan ( 3 )  -  tan (a) 
D*tan(a)tan(8) 


(2) 


y  = 


tan(3)  -  tan(a) 

Since  the  reported  angles  of  bearing  are  random  variables,  and  each 
pair  of  reported  angles  uniquely  describes  a  point  in  the  x,y  plane,  we  can 
transform  the  density  from  the  o,8  plane  to  a  density  in  the  x,y  plane, 
assume  the  joint  density  of  A  and  3  is  given  by  equation  (1),  f(a,B). 
define  two  new  random  variables  X  and  Y  as  functions  of  A  and  B: 

D  *  tan ( 3 ) 


X  = 


Y  = 


tan(B)  -  tan(A) 
0'tan(A)tan(3) 


(3) 


tan(B)  -  tan ( A ) 

Equivalently,  A  and  B  can  be  written  as  functions  of  X  and  Y 
A  =  Arctan[Y/X] 

3  =  Arctan[Y/ (X-0) ] 

The  joint  density  of  X  and  Y  is  then  simply 

g ( x , y )  =  f{Arctan[y/x] ,Arctan[y/(x-D) ]}• |J| 
where  |J|  is  the  absolute  value  of  the  Jacobian 


(4) 


(  =  ) 


J  - 


6a  6a 
6x  5y 
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6x  6y 


(5) 


The  density  is  then 


g(x,y)  = 


2tto 


A°3 


y-D 


(xl+yl )  [  (x-D) :  *-y*  ] 


(?) 


C  [Arctan(y/x)-Arctan(y  /x  )]*  [Arc tan ( y/ ( x-0) ) -Arctan ( y  / ( x  -D) ) ] ] 

expj^-j  ( -  +  - — - - - 
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where  ( x 0 , y0 )  is  the  intersection  of  the  lines  of  bearing.  It  corresponds 
to  the  point  («0,B0) . 

Undoubtedly,  this  density  function  makes  meaningful  analysis  difficult. 
Therefore  some  other  attack  is  called  for.  In  the  following  sections,  we 
will  describe  and  compare  different  approaches  to  the  two-sensor  problem. 
This  will  include  the  classical  metnod  devised  by  R.G.  Stansfield,  and 
numerical  techniques  devised  by  the  clinic. 
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III  CLASSICAL  APPROACH 

Last  year's  clinic  reported  tnat  the  classical  approach  to  finding 
probability  regions  yielded  ellipses  which  contained  only  approximately  50% 
to  80%  of  the  probability  claimed.  Because  of  the  magnitude  of  this 
discrepancy,  we  checked  this  method  to  see  if  the  previous  clinic  had  made 
an  error.  We  also  needed  to  understand  the  classical  approach  so  that  we 
could  compare  it  to  other  methods  of  analysis. 

The  classical  approach  was  developed  primarily  by  R.G.  Stansfield 
during  World  War  II  and  published  in  1947  (reference  1).  Because  of  the 
complexity  of  the  density  function  (7),  and  because  he  obviously  did  not 
have  access  to  computers,  Stansfield  made  a  simplifying  assumption  of 
parallel  displacement  of  lines  of  bearing.  This  assumption  translates 
angular  error  to  distance  error  and  allows  the  problem  to  be  modelled  using 
a  bivariate  normal  distribution. 

His  measure  of  error  is  the  distance  from  the  true  location  of  the 

* 

emitter  (or  target)  to  the  reported  line  of  bearing.  It  is  defined  as  the 
length  of  a  line  segment  per  pend i cu 1 ar  to  the  crue  line  of  bearing  that 
meets  tha  reported  line  of  bearing  (see  figure  3).  This  distance  is  labeled 
PI.  The  distance  from  the  true  location  of  the  emitter  to  the  sensor  is  DI , 
and  the  angular  error  is  P . 

The  standard  deviation  in  Stansfield's  work  is  estimated  using  the  fact 
that  for  small  angles  sin(v)  is  approximately  equal  to  P.  Also,  the 
distance  01  can  be  estimated  with  reasonable  accuracy  by  the  distance  01 ' 
(see  figure  3).  By  the  simple  trigonometry  of  the  model,  we  k.iow  that 

PI  *  01  x  sin(iji)  (8) 

which  can  be  estimated  as 

PI  =  DI'  x  *  (9) 


r* 

(2 


awaJ 


Stansfield’s  standard  deviation  is  then  the  distance  01'  times  the  standard 
deviation  of  the  angular  reading  (in  radians). 

To  carry  this  model  to  the  two-sensor  case,  errors  PI  and  P2  will  be 
associated  with  each  of  the  two  sensors.  Stansfield  assumes  that  these 
distance  errors  are  independent  and  that  each  error  can  be  modeled  as  a 
normal  probability  distribution  about  a  mean  of  0  and  a  standard  deviation 
described  above.  Therefore,  we  have  a  joint  normal  probability  distribution 
with  the  random  variables  PI  and  P2.  The  density  function  is 

1  r  Pi  *  p2* 

f  (pi  ,p2)  3 - exp  | -i  ( —  ♦  — )|  (10) 

2ttoj02  l  02 1  J 

We  now  need  this  distribution  written  in  terms  of  points  rather  than 
the  distances  PI  and  P2.  Here,  Stansfield  makes  use  of  his  assumption  of 
parallel  displacement  of  the  lines  of  bearing.  With  the  distance 
measurements  PI  and  P2,  there  is  associated  a  unique  paral lelogram  with  one 
vertex  at  the  true  location  of  the  emitter  (see  figure  4).  For  each  pair  of 
distances,  there  is  also  a  unique  point  which  is  the  vertex  opposite  the 
true  location.  Hence,  the  distances  can  be  converted  into  points  with  (x,y) 
coordinates  (see  figure  4).  Again,  making  use  of  the  trigonometry  of  the 
model,  we  can  rewrite  PI  and  P2  in  terms  of  the  opposite  vertex  points 

(x.y). 


PI  3  -xsina  ♦  ycosa 
P2  3  -xsinfi  ♦  ycoss 


(11) 


With  this  conversion  formula,  Stansfield  transforms  the  bivariate 


normal  density,  replacing  PI  and  P2  and  multiplying  by  the  appropriate 
Jacobian  determinant.  The  exponential  term 


thus  becomes 


1  1 

- [(-sina)x  ♦  (cosa) y] 1  *  - [(-sinS)x  ♦  ( cos S ) y] 2 

V  °2* 

This  can  be  more  easily  handled  by  setting 

sin’a  sin’s 
\  =  -  ♦  - 

V  °2* 

sina  cosa  sins  cosS 

V  =  -  ♦  - 

V  °2* 

cos ’a  cos’s 


(12) 


vj  =  -  +  - 

°1*  o2’ 

so  that  the  exponential  term  can  now  be  written  as 
A x *  -  2 vxy  ♦  wy* 

The  Jacobian  determinant  may  also  be  simplified  and  written  as 
Xu  *  vJ 

Having  completed  the  transformation  from  distance  measures  to  points,  the 
joint  normal  density  now  looks  like 


f~\\i  *  v1 

g(x,y)  *  -  exp[-J(Xx*  -  2vxy  ♦  uy')]  (13) 

2it 

This  probability  density  can  be  put  into  a  more  convenient  form  by  a 
rotation  of  axes  to  eliminate  the  'xy1  term.  Setting 

■2v 

tan(20)  =  _  (14) 

\  -  u 


the  rotation  is  given  by 

x  =  x’cos0  -  y'sinfl 


(15) 


y  =  x'sInC  ♦  y'cos0 

Note:  where  Stansfield  uses  capital  X  and  Y,  we  are  using  x'  and  y*  to  avoid 
confusion  between  random  variables  (X  and  Y)  and  variables  denoting  a 


& 


rotation  of  axes. 


If  these  values  for  x  and  y  are  substituted  back  into  the  original 
equation,  the  result  is 

1  1 

Xx2  -  2vxy  +  uy2  =  _ x'2  +  _ y*  *  (16) 

a2  b2 


where 


1 

-  =  .\cos20  -  2yCos2  sinC  *  usin20 

a2 

1 

-  =  ,\sin23  ♦  2vcos3  sin0  *  ucos20 

b2 


(17) 


At  this  point  in  Stansfieid's  derivation,  an  error  appears  in  the  original 
publication.  Please  see  the  appendix  for  a  discussion  of  the  error. 

For  this  transformat  ion ,  which  was  unaffected  by  the  error,  the 
Jacobian  determinant  conveniently  reduces  to  1.  Thus,  after  two 
transformations,  one  to  switch  from  distances  to  points,  and  another  to 
rotate  axes,  we  wind  up  with  a  joint  normal  probability  distribution  with 
respect  to  x',y‘  coordinates,  whose  density  is  given  by 


-  V2  r  1  1  -j 

h(x',y1)  =  -  exp  -}( - x'2  ♦  - y’2)  (13) 

2  it  L  a2  b2  j 

Stansfieid's  assumption  of  parallel  displacement  allows  him  to  make 
these  transformations.  Without  it,  there  would  be  no  meaningful  way  to  move 
from  the  distance  errors  PI  and  P2  to  x,y  coordinates.  The  resulting 
expression  is  the  relatively  tractable  bivariate  normal  distribution.  In 
fact,  one  of  the  properties  of  the  bivariate  normal  is  that  the  most 
efficient  probability  regions  will  be  perfect  ellipses  determined  by  the 
equation  (next  page) 
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1  1 

- x'1  ♦  _ y '*  =  -2ln ( 1 - p )  (19) 

a1  b’ 

where  p  is  the  level  of  probability  desired. 

As  noted  in  the  appendix,  the  formulas  for  the  coefficients  a  and  b 
above  are  printed  incorrectly  in  Stansfield's  paper.  The  incorrect  version 
makes  the  values  of  a  and  b  larger  than  they  should  be.  Thus,  ellipses 
generated  from  the  incorrect  version  would  be  too  small  and  capture  far  less 
probability  than  claimed.  This  would  explain  the  discrepancy  in  last  year's 
results . 
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IV  CLINIC  ANALYSIS  OF  STANSFIELO  ELLIPSES 

With  Stansf  ield '  s  density  (18),  the  problem  of  finding  a  probability 
region  is  greatly  simplified.  For  any  two  angle  readings  a  and  8,  and  a 
desired  level  of  probability,  p,  an  ellipse  can  be  formed  which  supposedly 
contains  the  point  of  intersection  of  the  reported  lines  of  bearing  with 
probabil ity  p.  However,  because  of  the  simplifying  assumption  of  parallel 
displacement,  the  ellipse  will  not  contain  exactly  that  amount  of 
probability.  The  clinic  therefore  tested  the  accuracy  of  the  Stansfield 
ellipse  as  a  probability  region. 

Given  a  region,  R,  in  the  x,y  plane,  and  given  a  two  dimensional 

density,  g(x,y),  one  can  determine  how  much  probability  is  contained  in  R  by 

evaluating  the  integral 
► 

g(x,y)  dA 

V 

The  clinic  performed  the  integration  letting  R  be  a  Stansfield  ellipse,  and 
g(x,y)  the  density  given  by  (7).  However,  since  we  could  not  analytically 
integrate  g(x,y),  we  used  numerical  methods  on  the  computer. 

To  verify  the  accuracy  of  the  numerical  results,  the  clinic  evaluated 
the  double  integral  using  both  an  existing  International  Mathematics  and 
Statistics  Library  (IMSL)  subroutine  (D8LIN),  and  a  modified  Simpson's 
integration  algorithm  written  by  the  clinic.  The  two  methods  yielded 
virtually  the  same  results. 

By  varying  the  levels  of  probability  over  three  sets  of  angle  readings 
with  o  *  2.0*,  we  consistently  found  tnat  the  amount  of  probability  claimed 
by  Stansfield  was  very  close  to  the  amount  actually  contained  in  an  ellipse 
(see  table  1  at  the  end  of  this  section).  Any  d  i  sc repenc  ies  were 
insignificant  and  may  be  partially  due  to  numerical  approximation  error. 
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However,  this  conclusion  holds  only  if  Stansfield  ellipses  as  corrected 
by  the  clinic  are  used.  Ellipses  derived  directly  from  Stansfield  '  s  paper 
(reference  1,  p.  763,  equations  (15),  (16),  and  (17)),  contain  only  about 
50?.  to  80S  of  the  probability  claimed,  which  is  what  last  year’s  clinic 
concluded . 

We  then  examined  the  difference  between  the  probability  in  the  top  and 
bottom  halves  of  an  ellipse.  Under  Stansfield ' s  assumption  of  parallel 
displacement,  the  probabilities  in  the  top  and  bottom  halves  are  equal. 
However,  upon  dropping  the  parallel  displacement  assumption,  we  found  that 
the  probabilities  contained  in  the  top  and  bottom  halves  are  actually  not 
equal.  Instead,  the  top  half  consistently  captures  less  probability  than 
the  bottom  half,  although  surpr isingl y,  the  sum  of  the  two  halves  come  very 
close  to  the  total  probabil ity  desired  (see  table  2). 

At  higher  probability  levels,  the  percentage  difference  between  tne  top 
and  bottom  halves  is  smaller  relative  to  that  of  lower  levels.  To  explain 
this,  we  examined  a  contour  map  of  the  clinic  density,  as  well  as 
probabil ity  bands  for  the  Stansfield  ellipse.  The  contour  map  (see  figure 
5)  graphically  demonstrates  that  the  probability  in  the  lower  half  is 
concentrated  nearer  the  center  of  the  ellipse,  while  the  probability  in  the 
upper  half  is  more  dispersed.  As  we  increase  the  desired  probability  level, 
the  rate  at  which  the  upper  half  of  the  ellipse  captures  probability 
increases  relative  to  the  rate  of  the  lower  halt. 

We  also  explored  this  issue  of  denser  bottom  half  probabilities  by 
looking  at  ellipses  in  terms  of  eight  probability  bands  of  equal  width  (see 
figure  6).  We  consistently  found  bands  close  to  the  center  below  the 
location  of  the  emitter  to  contain  more  probability  than  corresponding  bands 
above,  while  the  opposite  was  true  for  corresponding  bands  further  away  from 


the  center  (see  table  3).  This  indicates  that  when  the  Stansfield  ellipse 
is  made  large  enough,  most  of  the  additional  probability  must  come  from  the 
upper  half. 

Although  Stansfield's  parallel  displacement  approximation  is  not  a 
strict  depiction  of  reality,  his  method  is  a  simple  way  to  find  a  region 
which  contains  a  desired  level  of  probability.  This  does  not  necessarily 
mean,  however,  tiat  his  regions  are  the  "most  efficient",  in  the  sense  that 
they  capture  the  desired  level  of  probability  in  the  smallest  possible  area. 
The  clinic  wanted  to  find  out  if  it  was  possible  to  capture  the  same  amount 
of  prcbabil ity  in  smaller  regions.  Sections  V  and  VI  describe  regions  we 
examined . 
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TABLE  1 

Angle  readings  are  given  as  o/3 

80/110  50/60  40/150 

Claimed  Actual  Actual  Actual 


Prob. 

Prob. 

%  Diff 

Prob. 

5  Diff 

Prob. 

%  Diff 

0.25 

0.2512 

0.46 

0.2622 

4.39 

0.2500 

- 

0.50 

0.5029 

0.53 

0.5287 

5.73 

0.5000 

- 

0.75 

0.75C4 

0.06 

0.7452 

-0.65 

0.7495 

-0.07 

0.95 

0.9396 

-1.09 

0.3791 

-7.46 

0.9490 

-0.10 

0.99 

0.9802 

-0.99 

0.9225 

-6.31 

0.9897 

-0.07 

Note:  o 

=  2.0*  and 

the  distance  between 

sensors 

is  1000. 

I 
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TABLE  2 


80/110 


Claimed 

Prob. 

Top 

8ottom 

Total 

%  Oiff 

0.25 

0.1162 

0.1350 

0.2512 

16.24 

0.50 

0.2279 

0.2750 

0.5029 

20.70 

0.75 

0.3419 

0.4085 

0.7504 

19.47 

0.95 

0.4475 

0.4921 

0.93S6 

9.96 

0.99 

0.4738 

0.5014 

0.9802 

4.73 

50/60 


Claimed 

Prob. 

Top 

Bottom 

Total 

%  Oiff 

0.25 

0.1006 

0.1616 

0.2622 

60.57 

0.50 

0.1903 

0.3384 

0.5287 

77.80 

0.75 

0.2828 

0.4624 

0.7452 

63.54 

0.95 

0.3811 

0.4980 

0.3791 

30.69 

0.99 

0.4220 

0.5005 

0.9225 

13.60 

40/150 


Claimed 

Prob. 

Top 

Bottom 

Total 

%  Oiff 

0.25 

0.1244 

0.1256 

0.2500 

0.94 

0.50 

0.2480 

0.2520 

0.5000 

1.62 

0.75 

0.3701 

0.3794 

0.7495 

2.54 

0.95 

0.4654 

0.4836 

0.9490 

3.91 

0.99 

0.4838 

0.5055 

0.9893 

4.49 

Note:  a 

=  2.0*  and 

the  distance 

between  the  two 

sensors 

TABLE  3 


80/110 

50/60 

40/150 

Top  Band 

0.01276 

0.01855 

0.00982 

2nd  Band 

0.05649 

0.05355 

0.05157 

3rd  Band 

0.15208 

0.11777 

0.15229 

4th  Band 

0.26707 

0.23219 

0.27010 

5th  Band 

0.28376 

0.33246 

0.27993 

6th  Band 

0.16456 

0.16390 

0.16453 

7tn  Band 

0.04573 

0.00416 

0.05293 

Bottom  Band 

0.00491 

0.00000 

0.00807 

Total  Prob. 

0.98742 

0.92258 

0.98929 

■Note:  o  *  1.0*  for  the  80/110  case  and  2.0*  for  the  other  two  cases 


V  CLINIC  REGION 


The  clinic  attempted  to  find  a  probability  region  that  was  more 
efficient  than  the  Stansfield  ellipse.  That  is,  we  wanted  to  pack  as  much 
probability  into  the  smallest  area.  We  hoped  to  be  able  to  do  this  when  we 
dropped  the  parallel  displacement  assumption.  But  as  we  have  seen,  when  we 
drop  this  assumption,  the  resulting  probability  density  function  is  very 
complicated.  Hence,  we  attacked  the  problem  using  a  different  method. 

We  first  examined  the  probability  of  a  sensor  returning  a  reading 
resulting  in  a  line  of  bearing  inside  a  specified  sector  or  "fan".  In 
figure  7,  the  true  bearing  is  a,  but  we  assume  the  reading  sensor  I  returns 
is  a  random  variable.  A,  distributed  normally  with  mean  a  and  variance 
Hence,  we  can  easi  1  y  determine  P(alsAsa2).  It  is  simply  P(Asa2)-P(A$al )  . 
And  since  we  assume  the  underlying  probability  distribution  is  the  normal, 
we  can  compute  close  approximations  to  these  probabilities  on  the  computer. 

Next,  we  did  the  same  for  the  second  sensor,  so  that  we  have  the 
probability  that  sensor  2  returns  a  reading  inside  the  fan  between  bl  and 
b2  (see  figure  3),  where  8  is  distributed  normally  with  mean  0  and  variance 

V- 

The  two  fans  intersect  and  form  a  quadri 1 ateral .  The  probability  that 
the  two  sensors  return  readings  which  intersect  inside  the  quadri 1 ateral  is 
simply  the  probability  that  sensor  1  returns  a  reading  inside  fan  1  and 
sensor  2  returns  a  reading  inside  fan  2.  We  assume  the  two  sensors  are 
independent,  so  to  find  the  probability  of  the  quadrilateral  ,  we  multiply 
the  probabilities  of  the  fans. 

The  next  step  was  to  form  the  3-o  fans  (see  figure  9).  That  is,  we 
found  the  fan  formed  by  moving  £3o  degrees  from  a,  and  the  fan  formed  by 
moving  i3o  degrees  from  8.  The  reason  we  use  these  3-o  fans  is  simple.  The 
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quadrilateral  formed  by  the  intersection  of  these  fans  contains  99. 5%  of  the 
total  probability.  To  capture  100%,  we  would  have  to  examine  the  entire 
plane,  which  is  simply  not  possible.  So  we  only  look  at  a  subset  of  the 
plane,  but  we  choose  a  subset  such  that  we  lose  only  0.5*.  of  the 
probabi 1 ity. 

We  broke  this  quadr  i  1  ateral  up  into  smaller  subquadri  1  ateral  s , 
recording  the  area  and  probability  of  each.  As  can  be  seen  in  figure  10, 
each  subquadr i  1  ateral  has  a  different  area.  That  means  that  while  two 
subquadr i  1 ateral s  may  have  the  same  probability,  their  densities  are 
different.  In  other  words,  the  amount  of  probability  they  pack  into  a  unit 
area  is  different. 

We  then  made  a  very  fine  grid  of  3600  subquadrilaterals  and  took  those 
of  highest  density  to  form  a  probability  region.  Since  they  were  of  highest 
density,  they  contained  tne  most  probability  in  the  smallest  area.  We  made 
a  list  of  tne  subquadr i 1 ateral s ,  ordering  them  from  highest  density  to 
lowest.  We  formed  a  probability  region  by  connecting  subquadr i 1 ateral s 
together,  starting  at  the  top  of  the  list  and  moving  down  until  we  reached 
the  desired  level  of  probability.  The  areas  of  these  subquadrilaterals  were 
added  together  to  determine  tne  total  area  cf  our  probability  region.  Since 
every  one  of  the  subquadrilaterals  in  tne  region  was  of  higher  density  than 
the  subquadri 1 ateral s  outside,  we  concluded  that  this  was  the  most  efficient 
probabil ity  region. 

In  plotting  our  region,  rather  than  finding  and  drawing  possibly 
hundreds  of  quadrilaterals  linked  together,  we  assigned  the  probabil ity  of  a 
subquadr i  1  ateral  to  its  center  point.  In  this  way  we  obtained  a  list  of 
points  to  represent  the  subquadr i 1 ateral s ,  The  probability  region  was  then 
the  convex  hull  of  these  points. 
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We  defined  the  center  point  of  a  subquadrilateral  as  the  point  of 


VI  IMAGE  REGION 


Another  approach  the  clinic  examined  was  a  mapping  of  the  most 
efficient  probability  region  in  the  a,s  plane  to  the  x,y  plane.  Our  results 
showed  that  this  approach  is  not  as  efficient  as  the  Stansfield  and  Clinic 
methods. 

Recall  that  we  have  assumed  that  the  readings  the  two  sensors  return 
are  random  variables,  A  and  3,  each  distributed  normally.  The  joint  density 
of  A  and  3  is  therefore 
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(21).  This  created  a  region  in  the  x,y  plane  that  resembles  an  ellipse. 
The  area  of  this  region  was  calculated  by  the  vector  formula  for  the  area  of 
a  polygon: 

Area  =  l(xiy2n2y3+. .  •  •  -ynxl )  (22) 

Since  we  are  approximating  the  area  inside  a  smooth  curve  by 
calculating  the  area  of  an  inscribed  polygon,  there  will  be  some  error  in 
the  calculations,  but  the  following  analysis  shows  that  this  error  is  small. 

We  examined  the  accuracy  of  the  vector  formula  (22)  in  approximating 
the  area  of  an  elliptical  region  by  generating  36  unevenly-spaced  points  on 
an  ellipse.  We  then  compared  the  area  of  the  polygon  formed  by  the  convex 
hull  of  these  points  using  (22)  to  the  area  of  a  true  ellipse,  which  is 
given  by  the  formula  »ab  where  a  and  b  are  the  lengths  of  the  axes  of  the 
ellipse.  These  points  were  obtained  by  mapping  equally-spaced  points  on  a 
circle  to  an  ellipse  using  Stansfield's  transformation  (see  reference  1, 
equations  9  and  10). 

When  we  used  36  points,  the  relative  error  of  the  approx i mat i on  was 
0.005.  For  360  points,  the  relative  error  was  0.00005.  Thus,  the  vector 
formula  for  the  area  of  a  polygon  gives  a  sufficiently  accurate 
approximation  to  the  true  area  of  the  Image  region. 
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VII  RESULTS 

The  following  tables  give  comparisons  of  the  Clinic  and  Image  regions 
with  Stansfield  ellipses.  In  the  instances  where  the  Clinic  region  did  not 
contain  exactly  the  same  amount  of  probability  as  the  Stansfield  ellipse, 
figures  for  regions  capturing  slightly  more  are  included.  The  Image  region 
always  contains  the  same  amount  of  probability  as  the  Stansfield  ellipse. 

We  determined  probabil ity  regions  ranging  from  10%  to  95%  for  each  of 
three  sets  of  c ./fl  readings,  80/110,  50/60,  and  85  /95  ,  for  sensors  located 
1000  units  apart.  In  all  cases,  the  standard  deviation  is  l*=it/130  Rad. 

for  the  tables,  the  columns  are  labeled  as  follows 
A:  Probabil ity  captured  in  the  Stansfield  ellipse 
B:  Probability  captured  in  the  Clinic  region 
C:  Area  of  the  Image  region 
D:  Area  of  the  Stansfield  ellipse 
E:  Area  of  the  Clinic  region 

F:  Percentage  difference  in  areas  of  Clinic  region  and  Stansfield 
ellipse  (1-0/E) 


TABLE  6 


85/95 

A 

8 

C 

0 

E 

F 

101 

0.101 

0.102 

38,684 

38,219 

35,987 

-0.062 

20% 

0.202 

0.202 

82,967 

80,946 

76,421 

-0.059 

30% 

0.304 

0.304 

134,344 

129,385 

123,621 

-0.047 

40% 

0.406 

0.407 

195,019 

185,304 

179,448 

-0.033 

50% 

0.508 

0.508 

263,154 

251,442 

246,176 

-0.021 

60% 

0.607 

0.608 

358,880 

332,388 

329,836 

-0.008 

70% 

0.704 

0.704 

476,015 

436,746 

436,252 

-0.001 

80% 

0.797 

0.797 

637,237 

583,830 

533,064 

-0.001 

90% 

.0.883 

0.883 

891,007 

835,272 

812,860 

-0.023 

95% 

0.926 

0.926 

1,114,062 

1,086,714 

1,014,349 

-0.071 

The  Image  region  was  always  larger  than  the  Stansfield  ellipse,  even 
though  it  contained  the  same  level  of  probability.  In  other  words,  the 
Image  region  was  never  as  efficient  as  the  Stansfield  ellipse.  Thus  we 
discontinued  our  analysis  of  this  region. 

On  the  other  hand,  we  see  that  the  Clinic  region  as  described  in 
section  V  is  virtually  always  smaller  than  the  Stansfield  ellipse,  even 
though  it  contains  the  same  level  of  probability  or  more.  Hence,  we  say  the 
clinic  region  is  more  efficient. 

However,  it  is  only  slightly  more  efficient.  The  Stansfield  ellipse  is 
a  simple  analytic  region  easy  to  compute  and  draw.  The  Clinic  region  is 
much  more  difficult  and  time-consuming  to  obtain  and  draw.  Since  the 
difference  is  so  small,  it  is  not  clear  whether  the  gain  in  efficiency  is 
worth  the  trouble. 

To  examine  this  issue,  we  looked  at  why  the  Clinic  region  was  more 
efficient.  We  feel  that  because  the  clinic  drops  the  parallel  displacement 
approximation  the  Clinic  region  builds  up  and  the  Stansfield  ellipse  builds 
out.  By  that  we  mean  the  Clinic  region  chooses  the  densest 
subquadrilaterals  to  fulfill  a  desired  level  of  probability.  Since  we  know 
that  the  probability  is  generally  denser  below  the  .rue  location  of  the 
emitter  than  above,  a  Clinic  region  contains  more  subquadrilaterals  from 
below  for  lower  levels  of  probability.  For  higher  levels,  it  continues  to 
append  subquadrilaterals  upward  after  most  of  the  lower  ones  have  been  used 
(see  figure  II).  Therfore  we  say  the  clinic  region  builds  "up". 

The  Stansfield  ellipse,  on  the  other  hand,  builds  "out".  For 
increasing  levels  of  probability,  it  simpl y  enc loses  the  lower  probability 
levels  with  larger  ellipses  and  performs  no  shifting  (see  figure  12) . 

As  a  result,  the  geometric  center  points  of  the  Clinic  region  move 


upward  for  increasing  levels  of  desired  probability  (see  figure  13).  In  an 
effort  to  make  the  Stansfield  ellipse  more  efficient,  we  artificially  forced 
It  to  build  up  rather  than  out.  We  did  that  by  relocating  its  center  to  the 
geometric  center  of  the  Clinic  region  at  the  particular  probability  level. 
We  then  compared  the  density  of  the  relocated  ellipse  to  the  density  of  the 
ellipse  when  it  was  centered  at  the  intersection  of  readings. 

The  results  were  insignificant.  The  Stansfield  ellipse  did  become  more 
efficient  but  the  gain  was  miniscule.  At  lower  levels  of  probability  we 
found  an  increase  in  efficiency,  but  at  higher  levels  the  densities  were 
nearly  equal.  However,  as  we  have  noted  earlier  in  Section  IV,  the 
probability  in  the  upper  and  lower  halves  of  tiie  Stansfield  ellipse  is 
always  unequall y  d i str ibuted . 

Thus,  while  the  Stansfield  ellipse  is  not  the  most  efficient  region,  it 
is  probably  very  close.  Also  since  it  is  so  simple  and  quick  to  find  and 
draw  the  Stansfield  ellipse,  it  is  a  very  useful  and  good  approximation. 
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VIII  THE  THREE-SENSOR  PROBLEM 

The  clinic  has  approached  the  three-sensor  problem  in  a  manner  similar 
to  the  two-?ensor  analysis.  Our  purpose  is  the  same,  we  want  to  find  out 
what  the  true  density  is  when  we  drop  Stansfield's  parallel  displacement 
assumption.  Our  method  of  analysis  was  to  discretize  the  continuous 
probability  density  by  generating  point  estimates  for  each  set  of  three 
angle  readings.  Then,  since  there  is  more  than  one  way  to  generate  point 
estimates,  we  examined  a  number  of  methods  and  used  them  to  obtain  discrete 
probability  densities  in  the  x,y  plane. 

In  the  two-sensor  case,  two  intersecting  fans  form  a  quadrilateral.  In 
the  three-sensor  case,  though,  three  intersecting  fans  will  virtually  never 
form  a  usable  subregion  (see  figure  14).  Therefore,  we  cannot  approximate 
the  continuous  density  as  readily  as  we  did  in  the  two-sensor  case.  We 
decided  to  discretize  the  density  as  described  in  the  next  paragraph. 

We  assign  all  tne  probability  of  each  of  three  fans  to  the  lines  of 
bearing  bisecting  them.  We  assume  the  sensors  are  independent,  so  the 
probabil ity  of  returning  a  particular  set  of  three  readings  is  the  product 
of  the  probabil ities  of  the  individual  readings.  The  three  lines  of  bearing 
form  a  triangle  (see  figure  15),  and  we  assign  all  the  probability  of  the 
set  of  readings  to  a  point  to  represent  the  triangle.  There  are  several 
ways  to  determine  a  reasonable  point  to  use,  and  we  examined  a  number  of 
possible  methods. 

Sta.nsfield  develops  a  maximum  likelihood  estimate,  hereafter  referred 
to  as  the  MLE,  as  a  point  estimate  within  a  triangle.  The  derivation  of  the 
MLE  proceeds  from  the  joint  normal  density,  which  Stansfield  uses  because  of 
his  parallel  displacement  assumption.  The  principle  of  maximum  likelihood 
estimation  states  that  for  a  given  set  of  observations,  in  this  case  the 
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angle  readings,  we  maximize  the  probability  density  witn  respect  to  the 
unknown  parameters  we  wish  to  estimate.  Here,  our  unknown  parameters  are 
are  the  x  and  y  coordinates  of  the  point  estimate. 

These  MLE  points  have  some  interesting  characteristics.  Within  the 
triangle,  they  are  generally  located  nearest  the  line  of  bearing  that  is 
closest  to  the  true  location  of  the  emitter,  and  they  are  farthest  from  the 
more  acute  angles  of  the  triangle.  Stansfield  mentions  this  in  his  paper, 
and  our  computer  analysis  appears  to  confirm  his  statement. 

We  also  looked  at  the  "center"  of  a  triangle  as  a  point  to  represent 
the  triangle.  There  are,  however,  a  number  of  ways  to  define  the  center  of 
a  triangle.  The  clinic  considered  five  different  approaches  to  determine  a 
point  to  represent  a  triangle  in  the  sense  that  the  probability  of  the  set 
of  readings  (determining  the  triangle)  is  assigned  to  that  point.  The 
approaches  or  methods  were  the  following. 

1)  Minimizing  the  sum  of  the  distances  from  each  vertex  to  the  center 

2)  Minimizing  the  sum  of  the  squares  of  the  distances  from  each  vertex 
to  the  center 

3)  Angle  bisection 

4)  Intersection  of  the  medians 

5)  Perpendicul ar  bisection. 

For  acute  triangles,  all  the  methods  gave  center  points  that  lay  inside 
the  region,  while  for  obtuse  triangles,  the  perpendicul ar  bisection  method 
left  us  with  a  center  point  located  outside  the  triangle.  Although  there  is 
only  a  one  in  four  chance  tnat  the  emitter  is  inside  a  triangle  formed  by 
any  three  lines  of  bearing  (see  reference  1),  the  triangle  is  nonetheless 
the  only  bounded  region  in  which  to  look  for  the  emitter.  Therefore,  since 
we  wanted  ths  representati ve  point  to  lie  inside  the  triangle,  we  did  not 
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further  consider  the  method  of  perpendicular  bisection.  Also,  it  minimizing 


sums  of  squares  (method  2)  and  intersecting  the  medians  (method  4)  gave 


identical  results.  Thus,  we  will  only  discuss  methods  1,  2,  and  3. 


The  first  method,  minimizing  the  sum  of  the  distances,  requires 


minimizing  the  function 


Dl(*,y)  a  /-(xrx)»  7  (yi-y)«  ♦  /-(x2-x)‘  *  (y2-y)1  -  f23) 

|/-(x3"x) *  ♦  (y3-y)* 

where  (x^ ,  y  3 ) ,  (x2,y2),  aid  ( X3 , y3 )  are  the  vertices  of  the  triangle.  Since 
minimizing  01  by  standard  calculus  methods  is  extremely  complicated,  a 


computer  program  was  written  which  searched  tne  triangle  for  the  x  and  y 


coordinates  for  which  the  function  01  was  minimized. 


The  second  method,  minimizing  the  sum  of  squares  of  the  distances, 


proved  to  oe  quite  simple.  It  required  minimizing  the  function 


D2(x,y)  =  (Xj-xJ'^yj-y)1  ♦  (*2~x)  *  «•(  y2-y) 2  ♦ 

(24) 

(x3_x),+(y  yy)* 

Taking  the  partial  derivatives  of  02  with  respect  to  x  and  y  and  setting 


equal  to  zero  yields  the  point  estimate 


(x,y)  =  [(xi*x2*X3)/3  ,  (yi+y2+y3)/3]  (25) 

The  last  method  is  angle  bisection.  It  required  finding  the 


intersection  of  the  three  lines  that  bisect  the  angles  of  the  triangle. 


Figure  16  shows  where  the  various  methods  put  the  representat  i  ve  point 


for  a  particular  triangle. 


The  clinic  tested  to  see  how  much  difference  the  selection  of  various 


representative  points  would  have  on  the  discretization  of  the  probability 


density. 


We  created  point  clouds,  associating  probabilities  with  each  point. 


These  point  clouds  were  created  by  varying  the  lines  of  bearing  from  each  of 
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the  three  sensors  by  Jo  between  *3o  and  -3o.  Thus,  for  each  sensor,  13 
lines  of  bearing  were  produced,  which  in  turn  created  13-^  or  2197  different 
triangles.  Then,  within  each  triangle,  we  calculated  the  representati ve 
points. 

For  purposes  of  comparison,  we  analyzed  two  candidates  for  the 
representative  point,  Stansfield's  MLE  and  the  geometric  center  defined  by 
the  sum  of  squares  of  the  distances.  We  used  the  sum  of  squares  method 
because  of  its  relative  ease  in  programming. 

In  order  to  determine  a  discrete  approximation  of  the  continuous 
density,  we  first  split  the  area  of  the  plane  containing  the  point  cloud 
into  a  40  by  40  grid  of  rectangles.  Within  each  of  the  resulting  1500 
rectangles,  we  tabulated  the  number  of  points  (if  any)  contained  and  sjmmed 
the  probabilities  associated  with  each  point.  This  sum  of  probabilities  was 
then  assigned  to  the  rectangle. 

To  calculate  the  size  of  a  region  at  a  particular  probability  level,  we 
summed  the  rectangles  of  highest  probability  until  we  reached  the  desired 
level.  The  area  of  the  probability  region  is  simply  the  sum  of  the  areas  of 
the  rectangles,  which  are  all  of  uniform  size. 

Our  results  show  that  the  choice  of  point  estimate  makes  a  great  deal 
of  difference  in  the  nature  of  resulting  point  clouds  and  in  the  areas  of 
corresponding  probability  regions  (see  figures  17  and  13  for  examples  of 
point  clouds).  The  area  comparisons  (tables  7  and  8)  show  this  clearly. 

Note  that  even  though  the  areas  obtained  using  the  MLE  points  are 
smaller,  it  does  not  necessarily  mean  that  the  method  provides  "best"  points 
to  choose  for  our  purposes.  It  is  not  certain  which  discretized  density 
most  accurately  reflects  the  true  density  in  the  x,y  plane. 

We  would  suggest  that  further  work  on  this  problem  be  directed  toward 
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finding  the  true  density  function,  as  was  done  in  the  two-sensor  case,  by 
means  of  a  transformation  that  would  map  the  three  angle  readings  tc 
coordinates  in  the  x,y  plane.  The  formulation  of  this  mapping  is  still 
uncertain,  since  these  different  point  clouds  represent  competing 
definitions  of  the  transformation  needed,  and  it  is  not  clear  which  of  these 
transformations  is  most  appropriate. 


T ABLE  7 

AREA  COMPARISONS 

Grid  region  created  from  MLE  points 
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1 
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Angle  readings 


Level  of  probabi 1 i ty 
0.9  0.75 


0.5 
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TABLE  8 

AREA  COMPARISONS 

Grid  region  created  from  "minimized  sum  of  squares"  points  (method  2) 


Angle  readings 


80,110,69 

50,60,41 

85,95,84 


Level  of  probabil ity 
0.9  0.75  0.5 


20630 

455157 

3668714 


11269 

273093 

1590094 


4334 

136546 

564417 


P 

i 


Angle  readings 


80,110,59 

50,60,41 

35,95,84 


Stansfield  el  1 ipses 

Level  of  probabil ity 
0.9  0.75  0.5 


80,110,69 

10877 

6594 

50,60,41 

220152 

132544 

35,95,84 

499097 

300430  1 

Percentage  difference 

Angle  readings 

Level 

0.9 

of  probabil ity 
0.75 

-47.3 

-51.6 

-86.4 


-41.5 

-51.5 

-32.2 


-24.5 

-51.5 

-73.4 


& 


In  these  cases,  the  Stansfield  ellipse  is  SMALLER  than  the  corresponding 
Grid  region  by  the  given  negative  percentages. 
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IX  SYSTEMATIC  ERROR 

The  previous  sections  of  the  report  discussed  random  errors  in  passive 
sensors.  This  section  examines  systematic  errors  in  active  sensors.  An 
active  sensor  reports  distance  as  well  as  direction. 

For  tnis  model  there  is  only  one  sensor,  which  returns  direction  and 
distance  readings  at  point  PI  (see  figure  19).  It  is  then  moved  a  distance 
0  to  point  P2,  and  returns  new  readings.  The  two  measured  direction  readings 
form  the  angles  a  and  8  and  the  two  distance  readings  originate  at  PI  and  P2 
and  are  labeled  01  and  02. 

We  assume  the  readings  the  sensor  returns  may  be  wrong.  As  an  example, 
in  figure  19,  it  is  clear  the  direction  and/or  distance  readings  are  off 
since  the  readings  at  point  PI  do  not  agree  with  the  readings  at  point  P2. 
We  assume  there  are  no  random  errors,  but  only  systmatic  errors  associated 
with  the  sensor  which  are  the  same  at  both  locations.  That  is,  if  the 
sensor  reports  a  direction  reading  which  actually  deviates  by  €  degrees  at 
point  PI,  it  will  report  a  direction  reading  deviating  by  €  degrees  at  point 
P2.  Likewise,  how  far  the  distance  reading  is  off  at  point  PI  ij  how  far 
tne  distance  reading  will  be  off  at  point  P2. 

We  first  examine  the  directional  error,  disregarding  distance  for  the 
moment.  In  figure  20,  we  see  that  if  a  and  t>  .re  the  true  bearings,  the 
emitter  will  be  located  at  point  E.  8ut  if  the  directional  error  is  €,  the 
emitter  will  be  located  at  point  E‘.  We  continue  in  this  fashion,  producing 
a  number  of  points  which  would  indicate  the  true  location  of  the  emitter  for 
several  directional  errors.  We  then  find  tne  Lagrange  polynomial  using 
these  generated  points.  We  now  have  the  equation  of  a  curve,  somewhere  on 
which  the  emitter  is  truly  located. 

Now  we  examine  distance  error.  In  figure  21,  we  see  that  if  D1  and  D2 
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are  indeed  the  true  distances,  then  the  emitter  must  be  located  at  either 
point  E  or  F.  But  if  the  distance  error  is  6,  the  emitter  would  be  located 
at  either  £'  or  F*.  We  continue  producing  points  which  would  indicate  the 
true  location  of  the  emitter  if  the  distance  error  were  various  given 
amounts.  We  assume  that  in  practice  some  points  can  be  disregarded,  such  as 
F  and  F',  because  it  may  be  clear  the  emitter  is  not  located  in  that 
direction.  We  then  find  the  Lagrange  polynomial  using  tne  generated  points. 
As  in  the  directional  error  case,  we  now  have  the  equation  of  a  curve, 
somewhere  on  whicn  the  emitter  is  truly  located. 

So  now  we  say  the  true  location  of  the  emitter  is  somewhere  on  curve  1 
(figure  22),  but  also  somewhere  on  curve  2.  Therefore,  the  true  location 
must  be  at  tne  intersection  of  these  two  curves. 
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APPENDIX 


The  CGS/CMC  Mathematics  Clinic  discovered  an  error  in  the  paper  by  R.G. 
Stansfield,  published  in  the  Journal  of  the  Institution  of  Electrical 
Engineering .  volume  94,  Part  IIIA,  1947,  under  the  title  "Statistical  Theory 
of  D.F.  Fixing". 

The  error  is  in  equations  (15)  and  (15),  p.  768,  which  give  the 
coefficients  for  the  elliptical  confidence  region  in  the  x,y  plane.  As  they 
stand  now,  the  equations  are 


(15)  _  *  2  \  -  'Jtanfl 


(16)  _  =  2u  +  vtan0 

b* 

They  should  be 


(15)  _  -  \  -  utan0 

a* 


(16)  _  =  v  *  vtanO 

b* 


The  mathematical  derivation  of  (15)  and  (16)  proceeds  as  follows 
Starting  with  (11),  p.  763, 


P ( x , y )  dxdy  =  —  /~\\ i  -  v*  exp[-}(\x*  -  2 vxy  ♦  uy1)]  dxdy 
2m 


Stansfield  transforms  the  density  to  a  system  of  X,Y  coordinates  rotated 
through  an  angle  0,  where  0  is  given  by 


tan  20 


The  transformation  is  given  by  (next  page) 
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x  =  Xcos0  -  Ysin0 
y  =  Xsin0  ♦  Ycos3 

1  1 

This  changes  \x*  -  2'Jxy  «■  uy*  into  _ X*  ♦  _ Y*  where 

a1  b* 

1 

-  =  xcos  * 0  -  2vcos0sin0  ♦  usin*0 

a* 

1 

-  =  xsin’3  ♦  2\>cos0sin0  *  ucos*0 

b' 

By  the  definition  of  tan  20,  we  know  that 
2v 

u  *  A  *  - 

tan  20 

1  2v 

so  that  -  =  \cos*0  ♦  (X  - )$inJ0  -  2vcos0sin0 

a*  tan  20 

2 

a  X(cos*0  *  sin*0)  +  vsin0(  - sin0  -  2cos0  ) 

tan  20 

1  -  tan?0 

=  x  ♦  vsin0(  - sin0  -  2cos3  ) 

tan  0 

2tan  0 

since  tan  20  =  _  .  Continuing 

1  -  tan!0 

1  COS0 

-  =  X  *  vsin0[(l-tan*0) { - )sin0  -  2cos0] 

a1  sin0 

=  X  +  vsin0[cos0(l  -  tanl0  -  2)] 

=  X  +  vsin0[-cos0( tan *2  *  1)] 

=  x  -  vsin0(cos0  $ec*0) 


=  X  -  vsin0(cos3 _ ) 

cos*0 


X  -  vtan0 
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The  computations  for  —  proceed  almost  identically,  except  we  replace 

b* 

x  instead  of  u 


X  =  u  - 


tan  20 


—  =  (u 
b* 


tan  20 


.)sin*0  ♦  ucos*0  ♦  2  cos3sin0 


=  u(cosJ0  +  sin20)  -  vsin0(  _ sin0  -  2cos0  ) 

tan  20 

1  -  tanJ0 

=  u  -  '->sin0(  _ sin0  -  2cos0  ) 

tan  0 


1  -  tanr0 

since  tan  20  -  _  .  Continuing 

2tan  0 


-  =  u  *  v si n0[ ( 1  -  tan*0)cos0  -  2cos0] 

b* 


=  u  -  vsin0  cos0(-tan*0  -  1) 
=  u  ♦  vsin0  cos0  sec*0 


=  u  ♦  vtan0 


Therefore,  we  assert  that  equations  (15)  and  (16)  should  be 


=  X  -  \»tan0 


(16)  -  =  u  +  vtan0 

b* 


+  ««•  *  «  «  «  *  m  ar«  m-«  to  «  w  «  Mummi 


fanl 


fan2 


5«wi  ssaiia 


ALPHA  -  50 
BETA  -  60c 


FIGURE  13 


CR  90Z  * 


X  CR  80Z 


•  ALL  STANSFIELD  ELLIPSES 
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X  CR  102 


A)  MINIMIZE  SUM  OF  DISTANCES 

B)  MINIMIZE  SUM  OF  SQUARED  DISTANCES 

(=  INTERSECTION  OF  THE  MEDIANS) 

C)  ANGLE  BISECTION 

D)  PERPENDICULAR  BISECTION 
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